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Abstract 


Direction  of  arrival  estimation  by  eigenstruetnre 
methods  requires  knowledge  of  the  array  covariance 
matrix  and  an  exact  characterization  of  the  array  in 
teems  of  gsometry,  sensor  gam  and  phase,  ate.  It  often 
happens  that  the  actual  sensor  gain  and  phase  are  per¬ 
turbed  from  their  assumed  nominal  values.  If  eigeaatruc- 
ture  methods  are  applied  with  incorrect  sensor  parame¬ 
ters,  the  method  essentially  breaks  down  or  at  beat  gives 
poor  results.  We  propose  a  new  approach  which  uses 
information  in  the  observed  covariance  matrix  to  correct 
for  these  effects.  This  method  yields  substantially 
improved  performance,  a  fact  illustrated  by  the  results  of 
computer  simulations. 


L  Introduction 

Eigenstrueture  methods  were  first  introduced  in  time 
series  analysis  for  extracting  harmonics  embedded  in 
white  noise.  They  were  generalized  in  a  significant  way 
by  Schmidt  (1080)  and  applied  to  the  direrVm  of  arrival 
(DOA)  estimation  problem.  Their  chief  advantages  over 
conventional  methods  are  that  they  produce  estimates 
possessing  apparently  higher  resolution  and  yield  asymp¬ 
totically  exact  estimates  for  the  DOA  and  other  signal 
and  noise  field  parameters.  These  methods  are  however 
known  to  be  sensitive  to  sampling  and  modeling  errors. 
We  are  concerned  here  with  modeling  errors  in  the  sensor 
gain  and  phase.  Sources  of  such  errors  include  the  sensor 
itself  and  its  related  electronics.  Sonar  arrays  are  a  typi¬ 
cal  example,  where,  errors  of  1  db  RMS  in  gain  and  6  * 
RMS  in  phase  may  be  routinely  expected. 

The  degradation  due  to  sensor  errors  has  previously 
received  some  attention.  Cox  (1073)  studied  them  in  the 
context  of  optimal  beamforming  and  several  others  such 
as  Quazi  (1982)  have  considered  their  effect  on  conven¬ 
tional  array  processing.  As  yet,  no  results  have  been 
reported  on  the  effect  of  such  errors  on  the  performance 
of  eigenstrueture  methods  for  the  DOA  estimation  prob¬ 
lem  nor  have  means  for  mitigating  these  been  proposed. 


In  this  paper,  we  propose  an  algorithm  for  estimat-  I 
ing  the  unknown  sensor  gain  and  phase  by  exploiting  the 
structure  of  the  array  covariance  matrix.  The  resulting 
estimatee  are  then  incorporated  into  the  array  model  to  i 
reduce  the  DOA  estimation  problem  to  the  standard  I 
form.  Simulation  results  are  presented  which  snbetan-  ; 
tkte  the  performance  of  our  method. 


EL  Problem  Formulation 

Like  all  problems  in  array  processing ,  we  can  formu¬ 
late  the  present  problem  in  a  wideband  or  a  narrowband 
version.  Here,  we  consider  only  the  narrowband  problem. 
In  the  following,  we  use  lower  and  upper  case  to  represent 
scalars,  lower  case  bold  font  to  represent  vectors  and 
upper  case  bold  font  to  represent  matrices. 

Consider  an  uniformly  spaced  sensor  line  array  with 
m  sensors.  The  uniformly  spaced  line  geometry  assump¬ 
tion  b  essential  to  our  method  for  the  phase  estimation 
problem  but  not  for  the  gain  estimation.  The  sensors 
have  an  unknown  gain  and  phase  response.  Assume  that 
there  are  4  narrowband  stationary  zero  mean  and  mutu¬ 
ally  uncorrelated  sources  centered  at  frequency  w0  and 
sufficiently  far  from  the  array  to  allow  a  planar  wavefront 
approximation.  Additive  noise  b  present  at  each  sensor 
and  b  assumed  to  be  stationary  zero  mean  random 
processes  and  independent  from  sensor  to  sensor. 

The  received  signal  of  the  i-th  sensor  r,  (l )  is  given 
by 


'•,(<)  =  v,f 


k  —i 


-I)  —  fia 


+  «,<<» 


where 


the  signal  emitted  by  the  k  -th  source 


f»  s=  the  direction  of  arrival  of  the  i-th  source. 


the  spacing  between  the  sensors 
the  speed  of  propagation 
the  additive  none  at  the  i -th  sensor 
the  gain  of  the  i  -th  sensor 


Appr'"--?  j  f 


■3  i  stri i.-  utioa  i-it  -  i . 

85  10  11 


14  1 


w5R»i 


Fi’f 


IVt'i 


4;  *  the  phase  of  the  •  -tk  sensor 

We  assume  4  <  m.  Rewriting  (1)  in  matrix  notation, 
we  have 

r(M-M{A.(<)  +  •(!)}  (Sa) 

where  r(t)  and  a(t)  are  the  m  X  1  vectors 

rr(0=*|r,(<)  •  rm(t)\, 

*r(0“l»i(«)  ••  ■„(*))  (2b) 

and  a  (I )  is  the  4  x  1  vector  of  impinging  wavefronts 
with 

•r(0=*|»  «(<)  •  »<(*)]  (2c) 

A  is  a  m  X  4  matrix  with  a  columns 
{•i  ,  k  »  1, ..  4 }  defined  by 

...  ie-r«d«-thj  (2d) 

and 

n  -7*»*  (2a) 

Therefore,  A  is  a  Vandermonde  matrix  whose 
are  the  steering  vectors  of  the  impinging  planar  wave* 
fronts.  Finally  V  and  ♦  are  m  X  m  diagonal  matrices 
with  the  diagonal  entries  given  by 

1*1-  -  4i  ,  1*1-  -  (2f) 

Multiplying  (2a)  by  its  conjugate  transpose  and  tak¬ 
ing  expectations  ,  assuming  the  additive  sensor  noises 
have  variance  o3,  and  are  uneorrelated  with  the  source 
signals,  we  can  write 

P  •  E  [  r  (I  )r  ♦{! )  | 

—  ♦♦{A  SAN-  e3 1  (3) 

where  P  is  the  m  x  m  covariance  matrix  of  the  array 
measurements.  We  note  that  since  ♦  and  ♦  are  diagonal 
matrices  with  real  and  complex  entries  respectively,  we 
have  the  relations 


**  -  *  .  **  ««"»/#(♦>, 

44*  =■  I  ,  ♦ 

(4) 

Defining  G 

»  we  have 

P  »  GRG*  - 

(5) 

where 

R  -  ASANrl 

(8) 

is  the  array  covariance  when  all  sensors  have  uniform 
gain  and  phase  and 

Q  -  ♦ASA*0*+o3I  =  *R  V  (7) 


■  the  array  covariance  when  all  sensors  have  uniform 
gain  but  with  possibly  varying  phase  4 ( . 

This  completes  the  definition  of  the  model. 

Our  approach  to  this  problem  is  predicated  on  first  : 
leaning  the  array  covariance  matrix  P  and  luca  extract* 


ing  from  this  the  sensor  gain  and  phase  {4i }  and  (4,  }■ 
We  then  use  these  estimates  to  find  the  covariance  (  as  it 
tuns  out  to  within  some  unknown  constants  )  for  a 
hypothetical  uniform  sensor  gain  and  phase  array. 
‘Hereafter,  the  normal  eigeastrueturo  methods  are 
directly  applicable  to  obtain  the  DO  As. 

m  Estimation  of  Sensor  Gaia  and  Phase 

From  observations  r  (t; )  at  time  ,  the  array 
covariance  matrix  is  estimated  by 

?  -  -Jr  EH*/  )*<«/)♦  (8) 

In  the  sequel,  in  the  interest  of  notational  convenience, 
we  use  P  instead  of  P  with  the  understanding  that  we 
have  only  an  estimate  for  P .  This  situation  will  motivate 
us  to  formulate  least  squares  solution  to  gain  and  phase 
determination  problems  rather  than  seek  exact  solutions. 
We  also  note  that  we  caa  make  P  arbitrarily  exact  by 
making  N  large. 

We  begin  with  the  sensor  pda  estimation  problem 
and  note  that  if  the  sensors  had  uniform  gain  and  phase,  , 
the  true  array  covariance  R  will  be  Toeplits.  This  fob  ' 
lows  from  the  assumptions  of  uniform  linear  array  and 
uneorrelated  sources. 

It  follows  that 

IIPl.,1  -  IHM**/  (®) 

with  |  -|  denoting  modulus.  Also,  since  we  have  only  an  i 
estimate  of  P ,  (0)  will  be  only  approximately  satisfied. 
Our  aim  here  is  to  exploit  the  Toeplits  property  of  R  to 
eliminate  |  [R 1  from  (0)  to  determine  i>i  ■ 

With  this  intent,  we  introduce  a  quantity  determined 
from  P  j 

F.j*  -  MIIPl.,1  /|(P  |*,|)  (10): 

If  we  choose  i,j,k  and  /  such  that  •-/  st-l,  Le., ! 

|  [P  ),,  |  and  |  (P  )m  |  lie  on  the  same  diagonal,  then  from 
the  Toeplits  property  of  R  and  (0),  (10)  we  get 

Mijl w  —  -ln*y  -  Ml  +1“*/  (II) 

Note  that  (11)  will  not  be  exactly  satisfied  due  to  errors 
in  estimating  P . 

Taking  all  such  (nonredundant)  relations  when  i,j  i 
and  k ,/  pain  that  lie  on  the  main/super  diagonals,  we  ] 

get  a  total  of  k+  =  £  {» ( *  -l)/2)  equations  of  the  type  • 
«— s 

(11).  We  can  write  these  compactly  as 

B  (lng>„  Intfrj,  ,  lnV>„  |r 

=  (••*.  #*,,«.  '  '  •  lr  .  •-/  =  b-f  (12) 

Where.  B  is  a  x  m  matrix  whose  rows  have  one  of  i 
the  following  forms:  | 

(a)  |0,  ..  ,0,2,0, ..  ,0,-2,0,  .. )  when  i  =■/  and  1»1.  All 

entries  in  this  row  are  sera  except  for  a  2  and  -2  at 
the  •  th  and  the  k  th  positions  respectively.  I 
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(b)  |0,  ..  ,0,1,0, ..  ,0,-l,0, ..  J  wken  «yff  «»d  >*=*  ' 
All  entries  in  this  row  are  xero  except  for  a  1  and  >1  j 
at  the  •  th  and  the  /  th  positions  respectively. 

(c)  l ..  ,0,1,0, ..  ,0,-l,0, ..  ,0,-1, 0, ..  ,0,1,0, ..  1  when 

< ,;  ,k  and  f  are  all  distinct.  All  entries  in  this  row 
are  xero  except  for  and  1  at  the  »th, 

jr'th.ir  th  and  the  /  th  positions  respectively. 

Examining  B  shows  that  the  1  X  m  row  vector 

(1,1,1. _ ,l|r  ia  the  lone  vector  spanning  its  row  noil 

space. 

We  can  solve  (12)  with  a  least  squares  approach 
using  the  singular  value  decomposition  of  B .  i.e., 

IMi.  W*  •••  (IS) 

—  BO |  •  •  •  •  •  •  |r 

Where  B  *  is  the  pseudo  invent  defined  in  terms  of  the 
n-1  singular  values  and  the  corresponding  left  and  right 
aingnlar  vectors  of  B .  Le., 

B*  —  i, »*,••» »Jr  (M) 

With  ■  ,■  and  ▼  ,•  being  the  left  and  right  aingnlar  veeton 
of  B  respectively.  E  ia  a  n-1  X  n-1  diagonal  matrix  of 
the  n-1  aingnlar  values  of  B .  We  note  that  (12) yidda 
the  minimum  norm  least  squarea  solution.  A  general  sola* 
tion  can  be  obtained  by  adding  an  arbitrary  aerial  times 
the  null  space  vector.  This  amounts  to  saying  that  we  can 
only  determine  the  sensor  grin  to  within  an  arbitrary 
multiplicative  constant.  j 

Now,  having  found  ,  we  can  use  (3)  to  determine 
Q .  Hence,  starting  with  the  measured  array  covariance 
matrix  P  in  which  both  gain  and  phase  were  unknown, 
we  have  found  the  matrix  Q  (  to  within  a  Malar  con* 
stant  )  -  a  covariance  matrix  with  uniform  sensor  gain  j 
but  with  the  phase  unknown. 

We  estimate  sensor  phase  by  noting  from  (7) 

angle  (|Q  ],y  )  =-  angle  ([R  ],y  )  +  d,  -  dy  ( 15) 

Where 

angle  (y )  =  aretan  (imag  (y  )/real(g )) 

Now,  defining 

*',/«  =  angle  ([Q  )  -  angle  ((Q  Jw  )  ( 18) 

Since  angle  ((Q  |„ )  is  xero,  (18)  is  meaningful  only  when 
i,j,k  and  f  lie  on  the  super  diagonals. 

Again  invoking  the  Toeplitx  property  of  R , 
angle  (|R  |,y )  depends  on  •  alone  and  it  follows  that 

«',Jw  =■  d,-dy-d»+df,  *-/=**-/  (171 

Taking  all  nonredundant  relations  from  the  super 
diagonals  excepting  the  singleton  (P  ],„ ,  we  get 

kg  *■  £  (i(i  -l)/2)  equations.  And  we  can  write  these 

i  -2 

compactly  as  i 


C  (dn  ds>  ds>  •  •  ■  •  ) r 

*  I  •  •  •  •  lr  .  *->  =  k'l  (») 

Where  [dt,  d*.  dx»  •  •  • .  d.)  is  a  m  Xl  vector  of  unk¬ 
nowns  (d<).  I  *  •  •  Mi#,  •  '  •  1  is  a  kg  X  1  vector  of 
aretan  differences  computed  from  Q  and  C  is  a 
i  kg  X  m  matrix  whose  rows  have  one  of  the  following 
I  forms: 

(a)  ( ..  ,0,1,0, ..  ,0,-2,0, ..  ,0,1,0, .. )  when  ;  —  k .  AH  j 

entries  are  xero  except  for  a  1,-2  and  1  at  the  | 
i  th,;  («d  )th  and  the  I  th  positions  respectively.  ' 

(b)  [ ..  ,0,1,0, ..  ,0,-l,0, ..  ,0,-l,0, ..  ,0,1,0, .. )  when 
i,j,k  and  I  are  all  distinct.  All  entries  are  xero 
except  for  a  1,-1, -1  and  1  at  the  ith,  ;th,kth  and 
the  fth  positions  respectively. 

Again,  by  inspecting  the  rows  of  C ,  it  is  dear  that 
there  are  two  row  vectors  vis  (l,l,l,...,l|r  and 
|l,2,3,...m|r  in  tbs  null  space  of  C .  The  first  null  spaee 
vector  implies  that  phase  can  only  be  determined  to 
within  an  arbitrary  reference,  an  matter  of  no  conse¬ 
quence  in  the  DOA  problem.  The  second  null  space  vec¬ 
tor  is  significant  and  implies  that  the  phase  can  only  be  i 
determined  to  within  an  arbitrary  progressive  phase  fac- 
1  tor.  This  in  turn  means  that  the  best  we  can  hope  to  do 
is  to  find  the  DOAs  to  within  an  arbitrary  rotation  fac¬ 
tor. 

Once  again,  we  need  to  find  the  least  squares  solu¬ 
tion  and  this  is 

Idt.  dj,  •  ‘  ,  d*lr 
-  C*  [  •••  m^,  •  lr  (10) 

With  C  *  being  the  pseudo  inverse  of  C  and  defined  as 
earlier  in  terms  of  the  singular  values  and  vectors. 

IV  Simulation  Results  1 

Computer  simulations  were  carried  out  to  verify  the 
performance  of  the  proposed  approach.  The  simulation 
model  consisted  of  an  8  element  array  with  half 
wavelength  inter  element  spacing.  Three  equal  power  sig¬ 
nal  sources  were  generated  to  arrive  at  angles  -35  \  -83  * 
and  -25  '.  Further,  additive  and  uncorrelated  sensor  noise 
was  injected  with  a  sensor  level  SNR  of  0  db  referenced 
to  the  signal  sources.  The  sample  covariance  P  was 
accumulated  from  300  snapshots  of  array  data.  Fig.  1, 
plot  A  shows  a  ’spectral'  plot  obtained  by  eigenstructure 
methods  for  the  known  (identical)  sensor  grin  and  phase 
model.  As  expected,  we  get  peaks  in  the  true  directions 
and  all  the  sources  are  clearly  resolved.  Plot  B  is  for  the 
case  when  sensor  gain  and  phase  are  perturbed  from  their 
assumed  uniform  values.  O.Sdb  and  8  *  (both  RMS)  gam 
and  phase  errors  respectively  were  used.  An  examination 
shows  peaks  skewed  from  their  true  directions  and  a  reso¬ 
lution  performance  that  is  clearly  poor.  This  plot  is  only 
representative  of  the  degradation  to  be  expected,  since  it 
is  sensitive  to  the  actual  grin  and  phase  error  values 


4 


<u«d.  It  is  elenr  that  eigenstrncture  methods  are  badly 
degraded  is  several  practical  situations  when  sensor  data 
is  not  known  exactly.  > 

Fig.  2,  Plot  A  is  the  'spectral'  plot  a/ter  the  sensor  i 
phase  and  gain  have  bees  estimated  as  per  procedure 
detailed  in  section  ID  and  then  plugged  into  the  model 
before  the  eigenstrncture  algorithm  was  applied.  The 
new  procedure  has  obviously  succeeded  in  resolving  the 
three  sources.  However,  a  rotation  has  evidently  occurred 
due  to  the  null  space  vector  described  earlier.  Plot  B 
represents  the  corresponding  result  after  the  rotation  uoe- 
ert amity  has  been  resolved  by  means  of  an  additional 
input  such  as  the  phase  of  any  one  sensor  with  respect  to 
another.  The  plot  now  has  three  peaks  whieh  occur  in  the 
true  directions.  No  bias  is  evident.  It  is  clear  that  sensor 
phase  and  gain  estimation  procedures  have  succeeded  in 
restoring  the  resolution  and  accuracy  of  the  eigens tinc¬ 
ture  methods. 


In  this  paper  we  have  shown  that  sensor  pin  and 
phase  emeu  significantly  degrade  the  performance  of 

eigenstrncture  methods  in  the  DOA  estimation  problem. 
This  can  result  ,  in  poorer  performance  than  the  krw  reso- 
lotion  but  more  robust  conventional  angular  spectrum 
estimators.  We  then  showed  for  a  uniformly  spaced  Use 
array  and  uocorrelated  sources,  there  is  enough  informa¬ 
tion  in  the  covariance  matrix  for  estimating  the  sensor 
gain  and  phase.  These  estimates  can  then  be  used  to  fully 
restore  the  performance  of  the  cigenstructure  method. 
However,  since  the  gain  and  phase  can  only  be  deter¬ 
mined  to  within  certain  unknown  factors,  source  powers  i 
Can  only  be  estimated  to  within  an  arbitrary  scale  factor.  , 
More  significantly,  the  estimation  of  the  DOAs  can  be 
only  be  made  to  within  an  arbitrary  rotation  factor. 
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